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平衡 装载 约束 下 的 车 辆 路 径 问 题 研究 


雷 定 具 1f， 宋 文 杰 ， 张 英 贵 
(中 南大 学 交通 运输 工程 学 院 , 长 沙 410075) 


摘 要 : 针对 车 辆 三 维 装 载 约束 下 的 车 辆 路 径 问题 (3L-VRP) 进 行 研究 ， 引 进 车 辆 的 平衡 装载 约束 ， 综 合 考虑 传统 的 
先进 后 出 、 局 部 支撑 、 脆 弱 性 等 约束 , 构建 平衡 装载 约束 下 的 车 辆 路 径 问 题 (BL-VRP) 模 型 。 针 对 模型 中 的 平衡 约束 ， 
提出 一 种 接触 面积 的 装载 算法 。 在 此 基础 上 ， 构 建 以 回溯 遗传 算法 (B-GA) 为 骨架 的 多 阶段 算法 框架 ， 对 车 辆 路 径 优 
化 进行 求解 。 研 究 结 果 表 明 ， 多 阶段 算法 不 仅 在 解决 3L-VRP 上 好 于 目前 已 有 算法 ， 同 时 对 BL-VRP 表现 优秀 。 提 
出 的 多 阶段 算法 为 解决 BL-VRP 问题 提供 一 条 参考 思路 ， 但 在 时 效 性 上 需要 进一步 完善 。 

关键 词 : 物流 工程 ; 平衡 装载 ; 车 辆 路 径 优 化 ; 回溯 遗传 算法 ; 多 阶段 算法 
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Research on vehicle routing problem under balance loading constrain 


Lei Dingyou', Song Wenjie, Zhang Yinggui 
(School of Traffic & Transportation Engineering, Central South University, Changsha 410075, China) 


Abstract: This paper took vehicle routing problem (3L-VRP) with three-dimensional loading constraints as the research 
Pe 一 object，introducing equilibrium loading constraint of the vehicle, considering the traditional advanced backwards, local 
F support, and fragility constraints, constructing the vehicle routing problem (BL-VRP) model with balanced loading 
constraints. Considering the equilibrium constraints in the model, it proposed a loading algorithm of contact area. On this 
basis, the method used a multi-stage algorithm framework based on Backtracking Genetic Algorithms (B-GA) to solve 
vehicle routing optimization. The research results showed that the multi-stage algorithm was not only better than the 
existing algorithms in solving 3L-VRP, but also effective for BL-VRP. Although multi-stage algorithm providing a reference 
idea for solving the BL-VRP problem, it needs further improvement in timeliness. 

Key words: logistics engineering; balanced loading; vehicle routing optimization; backtracking genetic algorithms; 
multi-stage algorithms 


0 引言 模型 ， 设 计 多 阶段 算法 ， 构 建 了 路 径 优 化 策略 、 货 物 装载 顺 
序 生成 策略 以 及 单车 货物 装载 策略 三 个 阶段 ， 从 而 有 效 地 解 
平衡 装载 约束 下 的 车 辆 路 径 问 题 (BL-VRP) 是 三 维 装载 ” 决 该 问题 。 
约束 下 车 辆 路 径 问 题 (3L-VRP) 的 一 个 重要 扩展 。Gendreau 等 a 
人 由于 2006 年 首次 提出 3L-VRP， 在 研究 物流 的 配送 中 不 仅 1 ”BL-VRP 的 描述 及 模型 的 建立 
CS 仅 只 优化 车 辆 路 径 的 以 及 配送 的 需求 ， 同 时 考虑 货物 的 长 宽 1.1 问题 描述 
高 、 装 载 特性 以 及 “先进 后 出 ”原则 。BL-VRP 在 3L-VRP BL-VRP 问题 可 以 描述 为 : 给 定 配 送 网 络 G(V,E) , 其 中 : 
的 基础 上 ， 进 一 步 研究 货物 的 重量 、 重 心 对 车 辆 装载 稳定 性 V={0,1,2,.…n) 是 za+l1 个 顾客 ;五 是 连接 顾客 的 路 径 。0 是 配送 
的 影响 。BL-VRP 问题 在 现实 生活 中 普遍 存在 ， 比 如 家 用 电 ”中心 ;由 2… 帮 是 却 个 顾客 。 Cj=01…m 代表 顾客 i 到 jj 的 
器 、 家 具 以 及 建筑 材料 的 运输 ， 由 于 货物 的 重量 较 大 ， 不平 “路程 值 。 有 一 组 配送 货车 ， 车 厢 的 尺寸 分 别 为 工 、 到 、 卫 ， 
衡 的 装载 布局 将 严重 威胁 车 辆 行驶 的 安全 。 额定 载重 D ， 空 载重 量 ps。。 单 个 顾客 i 的 货物 总 数 m， 第 k 
前 国内 外 的 研究 结果 大 致 分 为 两 个 方面 ， 一 方面 对 ”个 货物 到 的 长 h、 宽 we、 高 hx、 体积 vx=lxxwxxhx 、 重 量 dx、 
3L-VRP 的 扩展 问题 进行 研究 ， 其 中 包括 时 间 窗 问题 P9、 车 底面 积 ax=lxwx 、 易 碎 性 所。 顾客 i 的 货物 总 体积 vol; ， 总 
回程 问题 申 、 同 时 取 送 货 问 题 仆 中等， 男 一 方面 的 对 基础 重量 weii 。 
3L-VRP 不 同 的 求解 算法 进行 研究 ， 有 精确 算法 办、 启发 式 算 把 车 厢 放 入 坐标 系 ， 车 厢 的 长 、 宽 、 高 平行 于 x 、》 
法 0 与 智能 优化 算法 。 由 于 3L-VRP 问题 由 两 个 NP- 难 问 ” 轴 。 货 物 到 在 最 后 最 左 最 下 的 顶点 坐标 为 Gu,ya,zx) ， 重 心 在 
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题 组 成 ， 解 决 难度 较 大 ， 求 解 算法 一 般 是 启发 式 算法 与 智能 货物 几何 中 心 。 

法 结合 。 其 中 禁忌 算法 022593 占 主要 地 位 ， 同 时 有 模糊 遗传 BL-VRP 约束 分 为 车 辆 路 径 约束 与 装载 约束 两 个 部 分 : 
算法 9、 蚁 群 算法 51 等。 现 有 研究 以 货物 尺寸 为 基础 ， 对 车 a) 车 辆 路 径 约 束 。 具 体 如 下 : 

辆 路 径 问 题 上 进行 扩展 , 较 少 涉及 货物 的 重量 、 重 心 的 影响 ， (R1) 货 车 在 配送 中 心 完 成 装载 ， 对 所 有 顾客 送 货 后 回 到 
无 法 使 货物 安全 有 效 地 配送 。 针 对 此 问题 , 本 文 提出 BL-VRP ”起 点 。 
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(R2) 每 个 顾客 只 由 一 辆 货 
(R3) 路 径 上 顾客 货物 总 重 
及 车 厢 容 积 。 
(R4) 每 条 路 径 的 装载 方案 必须 可 行 。 
b) 三 维 装载 约束 。 具 体 如 下 : 
(C1) 货 物 不 能 重合 ， 不 能 位 于 车 规 外 。 
(C2) 货 物 摆 放 时 ， 三 边 平行 于 车 厢 。 
(C3) 货 物 垂直 向 上 , 可 在 水 平面 旋转 。 假 定货 物 原始 长 、 


车 
量 、 总 体积 必须 小 于 货车 的 载 


加 高 分 别 为 Lo 、 Wio、 jako 。 


(C4) 先 进 后 出 约束 : 伙 载 时 ， 需 要 和 载 的 货物 不 能 被 之 


后 和 卸载 的 货物 阻挡 。 

(C5) 货 物 摆 放 时 ， 货 物 的 底部 需要 一 定 的 支撑 面 维 持 稳 
定 。 支 撑 面 由 车 厢 底 部 或 其 他 货物 的 顶部 构成 。 支 撑 面 要 达 
到 货物 底面 的 一 个 比例 。 


制 。 


衡 为 


(C6) 易 碎 货 物 上 面 只 能 摆 放 易 碎 货物 ， 非 易 碎 货物 无 限 


I 


(C7) 平 衡 约束 : 平衡 分 为 纵向 平衡 与 横向 平衡 。 纵 向 3 
轴 重 约束 ; 横向 平衡 研究 较 少 ， 本 文采 用 铁路 的 标准 ， 


化 
页 咎 


hl 


的 重心 横向 偏 移 中 线 面 不 能 超过 0.1 m。 


于 BL-VRP 的 约束 较 多 ， 满 足 上 述 约 束 的 路 径 解 ， 需 
要 得 到 相应 的 货物 装载 布局 方案 ， 从 而 确定 车 辆 的 稳定 性 ， 
同时 使 配送 行驶 的 总 路 程 值 最 小 。 
1.2 ”模型 的 构建 
于 模型 中 使 用 的 符号 较 多 ， 本 文 将 以 表 1 形式 给 出 。 
表 1 模型 参数 的 描述 
Table ] Description of model parameters 
符号 意义 符号 意义 
n 总 顾客 数 K 使 用 车 辆 数 
L 车 厢 长 度 (m) W 车 厢 宽 度 (m) 
末 车 厢 高 度 (m) D 额定 载重 (t) 
Do 空 载重 量 (t) Cy 顾客 i 到 jj 路 程 值 
mi ”顾客 i 的 货物 总 数 Ix 顾客 i 的 第 个 货物 
lix 货物 的 长 (m) Wik 货物 及 的 宽 (m) 
hix 货物 及 的 高 (m) Vik 货物 天 的 体积 m) 
dix 货物 及 的 重量 (t) Ci 货物 的 底面 积 (m?) 
fi 货物 的 脆弱 性 VOL 顾客 i 货物 总 体积 
weii ”顾客 i 货物 总 重量 Gxinyyin,zin) 货物 fx 的 左下 角 
N ”货车 /服务 顾客 数 a 货车 的 第 7 顾客 
根据 以 上 约束 , 以 货车 行驶 的 总 路 程 值 最 小 为 优化 目标 ， 


BL-VRP 数学 模型 下 : 
obj minf= > (Cow +D ,Cia n+ Cw0) (1) 
St. Xz Xr(rell,K] j,kell,N],jzk) (2) 
Xi Kli, j el[l, KJ,izj uell,Ni] vell,N)) (3) 
Nr>1(r ell,K]) (4) 
DN=n (5) 
vol:= > vali e[l,n]) (6) 
> volw < LxWxH(r ell, KD) (7) 
weii:= > dali e[l,n]) (8) 
bp weix; < D(r el1, K]) (9) 


Xik + li < XK Lit + hir < Dir Wir + Wik < 其， 


(10) 
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i,jel[ll,n] ke[ll,m] k'e[ll,m] hk 先 于 如 ' 放 
Xxtlx<LMzxthx<HAyrtwr<sW 


ie[l,N:] kell,mi] r ell,K] 


(lix, Wik, hit) = (liro, Wiro, hiro) J (wiro, liro, hiro) 
ie[l,n] kel[l,mi] 


Kik + lir < Kk Zit + hi < Tx Vir + Wik < Ye' 


i,j el[l,n] kell,mi] k'ell,m] 顾客 j 先 于 i 服务 


ax x_min max ymin y) > Odir 
(max xX=max(xxt', xix) min x= min(xx, xi)) 
(max y= max(yix, yx) min y= min( yt, yx)) 
(i,je[l,N:] rell,K] k'ell,m] kell,m] I A 处 于 Ix 
frzly fr=0, fr'=0 
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式 (10) 


(11) 


(12) 


(13) 


上方) (14) 


Gi,j elLN] rstK] kK'eltbm] kelbm] 到 :处 于 五 正 上 方 ) (15) 


区， di x (Xi + 3) > 


Re (16) 

ZLx> DY dr—NixL(r ell,K)) 
p23 dx Gt hy) < Nix Lr ell,K)) (17) 

EOet D+ DW wl 

5D a+ Do 2 
(N ell,N']) (18) 
标 函 数 式 (1) 求 总 路 程 的 最 小 值 ; 式 C2)~(4) 确 保 每 辆 货 
车 至 少 对 一 个 顾客 送 货 ， 每 个 顾客 只 送 货 一 次 ; 式 (5) 表 明 货 
车 服务 的 顾客 总 数 等 于 总 顾客 数 ， 式 (6)~(9) 确 保 顾客 货物 的 
总 质量 与 总 体积 小 于 货车 限制 ; 式 (10)(11) 阅 述 货物 之 间 不 能 
重 又 ,不 超过 车 厢 范 围 ; 式 (12) 表 明 货 物 可 在 水 平面 90°? 旋 转 ; 
式 (13) 保 证 满足 先进 后 出 原则 ; 式 (14) 表 示 稳 定性 约束 ， 选择 
支撑 比例 9 为 0.75; 式 (15) 描 述 货物 易 碎 性 约束 ; 式 (16)(17) 
为 轴 重 约束 ， 式 (18) 表 示 货 车 的 重心 在 装载 时 ， 偏 移 中 线 面 


的 值 不 超过 0.1 m。 
2 ”BL-VRP 问题 的 多 阶段 求解 算法 框架 


针对 BL-VRP 问题 ， 本 文 提出 多 阶段 优化 算 
基本 结构 如 图 1 所 示 。 该 算法 由 路 径 优化 决策 、 
序 生成 决策 、 单 车 货物 装载 决策 三 个 阶段 组 成 。 


随机 生成 种 
群 


初始 解 


| 坦 一 -| 
BGA 算 法 | | 人 
| 贪心 算法 


| 接触 面积 法 


性 入 顺 序 下 
成 决策 


| 


| 一 装载 决策 


图 1 


多 阶段 算法 结构 


Structural of multi-stage algorithms 


Fig. 1 


法 。 算 法 的 
货物 装载 顺 
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录用 定稿 
2.1 路 径 优化 决策 

随机 生成 的 初始 解 作 为 回溯 遗传 B-GA) 算 法 的 输入 ， 通 
过 单车 的 载重 和 体积 约束 ， 分 配 每 辆 货车 需要 服务 的 顾客 ， 


雷 定 就 ， 等 : 平衡 装载 约束 下 的 车 辆 路 径 问 题 研究 


ChinaXiv 合 作 期 刊 

第 37 卷 第 6 期 

小 依次 排列 。 
1) 接 触 面积 定义 

由 于 货物 可 在 水 平面 90° 转 动 ， 货 物 有 两 种 捍 放 方法 。 


形成 多 个 小 规模 TSP 问题 ; 然后 通过 回溯 法 遍历 ， 求 得 单车 
最 佳 行驶 路 径 以 及 路 程 值 。 
本 文中 采用 的 遗传 框架 与 思想 可 见 文献 [16]， 其 中 对 适 
应 度 的 计算 进行 说 明 。 本 文 以 货车 行驶 的 路 程 值 为 适应 度 ， 
每 辆 货车 的 路 程 用 回 济 算 法 计算 。 得 到 每 辆 货车 服务 的 顾客 
后 ， 回 渊 算法 寻找 该 货车 的 最 佳 行驶 路 径 以 及 路 程 值 。 回 漳 
法 是 一 种 深度 搜索 法 ， 按 照 优化 目标 进行 深层 搜索 。 当 搜寻 
到 某 一 步 时 ， 发 现 当前 的 选择 并 不 优 或 达 不 到 三 维 装载 检测 
要 求 , 退 一 步 重新 选择 。 具体 的 算法 思想 与 流程 参见 文献 [17]。 
2.2 货物 装载 顺序 生成 决策 
得 到 货车 行驶 路 径 后 ， 通 过 贪心 算法 ， 改 变 顾客 货物 的 
内 部 顺序 ， 形 成 不 同 的 货物 装载 序列 ， 调 用 装载 算法 。 如 果 
能 全 部 装载 ， 则 为 合法 路 径 ， 和 否则 路 径 非 法 。 
算法 1 贪心 算法 
a) 输入 货物 原始 的 装载 序列 了 7， 和 迭代 次 数 ier = 货物 个 数 ， 
nolmproved =0 ，i=1， totalVolume 是 货物 总 体积 ; 
b) 调用 算法 2， 得 装载 体积 packyoure 。 若 
packVolume=totalVolume ， 输 出 装载 成 功 ， 否 则 进入 步骤 c)。 
c) 每 个 顾客 货物 按 体 积 从 大 到 小 排列 , 形成 新 的 序列 。 
调用 算法 2， 新 的 装载 体积 packVolume' 。 
d) 若 
packedVolume = packedVolume' ， 


packedVolume' > packedVolume 则 
序列 7 六 代 蔡 。 若 
packedVolume=totalVolume ， 输 出 装载 成 功 ， 否 则 进入 步骤 e)。 

e) 改变 顾客 货物 的 内 部 顺序 产生 新 的 序列 二 。 再 次 调用 
算法 2， 得 装载 体积 packedVolumne' 。 
D 若 


packedVolume = packedVolume' , I 


packedVolume' > packedVolume 则 
7' 代替 , 重 置 nolmproved =0 。 


输出 装载 成 功 ， 否 则 


若 packedVolume=totalVolume ， 


nolImproved = nolImproved +1 。 


g) 返回 e)， 直 到 nolmproved =iter+1 。 

h) i=i+1， 重 复 上 步骤 b)~g)， 直到 i=5 。 若 未 输出 装载 
成 功 ， 则 输出 装载 失败 。 
2.3 单车 货物 装载 决策 

剩余 空间 指 货物 装载 后 剩余 的 未 使 用 空间 ， 列 表 
5S={s452.…5)} 表示 。 本 文采 用 最 大 空间 法 来 表示 剩余 空间 ， 
如 图 2 所 示 。 


图 2 最 大 空间 的 产生 
Fig.2 Generation of maximum space 
最 大 空间 指 由 货物 表面 或 车 厢 壁 与 车 厢 门 之 间 形 成 的 最 
大 矩形 空间 。 需 要 注意 的 是 ， 其 中 的 一 个 面 必 须 是 车 厢 门 。 
更 新 列表 5s 时 ， 列 表 按照 空间 坐标 的 了 最 小 、* 最 小 、z 最 


为 了 得 到 更 优 的 摆 放 方式 ， 本 文采 用 货物 与 空间 的 适应 度 
CS(b,s) 确定 待 装 货物 的 方向 以 及 相应 的 装载 空间 。5CS4w,9 表 
示 货 物 b 放 入 空间 ss 后， 与 周围 货物 或 者 车 厢 接 触 面 积 和 。 
Cs(b,s) 值 越 大 ， 货 物 与 空间 的 适应 度 越 高 ， 最 高 适应 度 的 空 
间 选 为 装载 空间 。 如 果 货物 不 能 放 入 空间 5 或 者 摆 放 后 支 
撑 底 面积 ai<0ax(96=0.75) ， 适 应 度 为 -oo 。 
2) 车 辆 平衡 约束 
本 文 对 车 辆 平衡 约束 分 为 纵向 轴 重 约束 与 横向 
约束 。 
a) 纵 向 轴 重 约束 。 
根据 车 厢 纵 向 受 力 ， 建 立 坐 标 系 ， 如 图 3 所 示 。 


Z 


考 
[5 
Em 
KS 


Q3=2N1 
CO4=NIi+N2 
图 3 纵向 重心 边界 条 件 分 析 


Fig.3 Boundary conditions analysis of longitudinal center of gravity 


Q@ 合 重心 位 置 在 (0,X) 时 ， 由 力矩 平衡 NixL=QOx(L-X1)， 


则 击 = 人 人 -X09 。 


~ 


@@ 合 重心 位 置 在 (XD 时 ， 由 力 和 


TOI 
| 


平衡 NxL=QOxX，， 则 


@ 合 重心 位 置 在 点 XX% 时 ， 装 载 最 大 重量 0=Ni+N2， 则 


b 模 向 重心 偏 移 约束 。 
横向 重心 偏 移 相对 于 纵向 加 
受 力 得 到 货车 合 


dx(Ow+wy)+TDoxW 


[wll 


约束 较为 简单 ， 根 据 横向 
心 的 横向 坐标 


[wl 
nh 


. pw di +Do 其 中 [WS|< 0.Im 即 可 。 
3) 接 触 面 积 算 法 
以 算法 1 产生 的 货物 序列 作为 输入 ， 输 出 可 装载 的 货物 
总 体积 。 具 体 流程 见 算法 2。 


算法 2 接触 面积 算法 

a) 输 入 货物 装载 序列 。 

b) 车 厢 的 体积 Va=LxWx 石 ， 初 始 列表 5= {Vn} 。 

Cc) 路 径 7 的 顾客 i=1， 顾 客 i 的 货物 k=1。 

gd) 遍历 列表 5 ， 计 算 货物 两 个 方向 的 适应 度 CS(5,s) 。 若 
CS(b,s)z-% ， 则 按 较 大 的 适应 度 值 对 应 的 货物 方向 放 入 相应 
装载 空间 最 后 最 左 最 下 的 角落 ， 更 新 空间 列表 


S={snspsp;， 车 C5(b,s)=-% ， 仓 载 已 装 的 顾客 ;的 货 


转 步 又 hh)。 
e) 对 列表 5 排序 ，k=k+1; 重复 步 又 d)， 直 到 k=m+1。 


作对 顾客 i 的 货物 进行 

中 ， 否 则 转 步 又 hh)。 
2)i=i+1， 重 复 步 骤 d)~fj， 直 到 i=N;+1。 
h) 输 出 已 装载 的 货物 总 体积 V 。 


3 ”实验 数据 及 分 析 


于 前 没有 BL-VRP 问题 的 相关 算 例 ， 本 文采 用 
3L-VRP 的 标准 算 例 以 及 物流 配送 实例 进行 两 次 实验 。 程 序 


平衡 约束 检测 ， 检 测 通过 ， 转 步 又 


以 Java 语言 进行 开发 ， 基 于 InteliJ IDEA 平台 运行 。 遗 传 算 
法 的 最 大 和 迭代 次 数 分 别 采 用 500、800、1 000、1 500， 种 群 


规模 采用 30、50、80、100 对 3L-VRP 标准 算 例 进行 运算 ， 
运算 时 交叉 概率 分 别 选取 0.7、0.8、0.9, 变异 概率 选取 0.05、 
0.08、0.1。 通 过 多 次 实验 ， 最 大 迭代 次 数 1 000、 种 群 规模 
50、 交 叉 概 率 0.8、 变 异 概率 0.05 的 结果 最 好 。 故 实验 中 各 
个 参数 设置 为 相应 的 值 。 
3.1 L-VRP 标准 算 例 实验 

使 用 多 阶段 算法 对 27 个 3L-VRP 的 标准 算 例 进行 10 次 
实验 ， 实 验 的 平均 结果 与 目前 的 已 有 的 VRLHI102 算 法 和 
TS-ILAD3] 算 法 进行 对 比 ， 如 表 2 所 示 。 表 中 dv 表示 多 阶段 

和 法 与 其 他 两 个 算法 最 好 结果 的 偏差 。 
表 2 3L-VRP 标准 算 例 实验 结果 
Table 2 Experimental results of 3L-VRP standard example 
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表 2 可 知 ， 与 YVRLHI1 算法 对 比 ， 每 个 案例 多 阶段 算 
法 得 到 的 结果 都 不 差 于 VRLHI1 算法 ; 与 TS-ILA 算法 结果 对 
比 , 在 顾客 数目 少 于 50 时 ,多 阶段 算法 稍 逊 于 TS-ILA 算法 。 
但 是 随 着 顾客 数目 的 增加 , 多 阶段 的 算法 结果 要 好 于 TS-ILA， 
算法 优越 性 得 以 体现 。 从 整体 上 来 看 ， 多 阶段 算法 的 运行 平 
均 结 果 要 好 于 其 他 两 种 算法 ， 说 明 本 文 提出 的 多 阶段 算法 在 
解决 3L-VRP 时 有 着 显著 的 效果 。 

3.2 BL-VRP 物流 配送 实例 

于 缺少 BL-VRP 的 标准 算 例 ， 本 文 对 长 沙市 某 地 域 配 
送 中 心 一 次 送 货 任务 进行 分 机， 并 进行 多 次 运算 取 最 好 值 ， 
同时 与 广泛 使 用 的 禁忌 搜索 算法 进行 比较 。 该 任务 有 15 个 顾 
客 ， 顾 客 位 置信 息 如 图 4 所 示 。 顾 客 之 间 以 及 顾客 与 配送 中 
心 之 间 的 距离 以 实际 的 最 短路 程 为 准 ,具体 信息 见 表 3( 单 位 : 
km) .顾客 的 货物 尺寸 信息 见 表 4, 其 中 每 个 顾客 货物 种 类 (Var) 


不 超过 三 种 , 每 种 货物 的 数量 (Num) 在 1~3 间 ( 长 度 单位 : m， 
重量 单位 ; t。) 


表 3 路 程 数 据 信 息 
Table 3 Route data information 
顾客 0 1 2 3 4 3 6 7 8 9 10 11 12 13 14 15 
0 0 le 1 8 125 17 1B I 9175 1 L221 233 21 3243 33 3716 
0 22.2 15.3 21.5 18.1 28.2 23.9 24.1 18.5 11.9 33.3 13.9 33.4.46.7 31.3 


要 0 16 20.1 24.7 10.4 10.7 24.1 31.5 31.1 15.6 29.6 23.8 21 50.1 
k 0 20.1 10.1 20.9 22.7 10.1 25.1 19.8 27.4 25.8 26.1 39.9 45.5 
4 0 29.5 33.5 8.8 27.9 10.3 37.4 28.8 18.2 33.5 22.8 27.4 
5 0 30.1 36.5 14.4 29.1 19.9 37 33.4 35 47.9 46.5 
6 0 18.429.434.837.9 9 36.5 17.5 43.2 55.9 
7 0 30.5 17.9 35 21.3 25.7 29.4 25.1 34.9 
8 0 34 27.3 34.4 38.4 39.7 47.7 51.4 
9 0 27.7 37.8 15.1 41.9 23.4 21.3 
10 0 44.2 16.6 43.3 55.4 26.3 
11 0 42.6 23.8 46.5 69.7 
12 0 44.2 39.7 16.6 
13 0 62.4 67.5 
14 0 45.1 
15 0 


算 例 名 称 VRLHI1 TS-ILA 多 阶段 算法 dv/% 
E016-03m 302.02 302.02 302.02 0.00 
E016-05m 334.96 334.96 334.96 0.00 
E021-04m 401.44 381.37 386.34 1.30 
E021-06m 437.54 437.19 437.19 0.00 
E022-04g 451.03 436.79 447.58 2.47 
E022-06m 498.38 498.32 501.06 0.55 
E023-03g 772.49 768.94 771.02 0.27 
E023-05s 821.35 805.77 808.55 0.35 
E026-08m 645.81 631.68 630.13 -0.25 
E030-03g 827.29 828.99 827.29 0.00 
E030-04s 815.62 780.61 778.05 -0.33 
E031-09h 630.46 614.6 610.23 -0.71 
E033-03n 2694.81 2636.85 2656.72 0.75 
E033-04g 1413.59 1398.77 1398.77 0.00 
E033-05s 1355.5 1352.76 1351.38 -0.10 
E036-11h 705.05 698.92 698.61 -0.04 
E041-14h 917.96 866.4 866.4 0.00 
E045-04h 1228.98 1228.47 1228.47 0.00 
E051-05e 753.87 763.09 750.17 -1.69 
E072-04f 596.42 590.99 79 -1.94 
E076-07s 1107 1096.53 1086.26 -0.94 
E076-08s 1171.49 1155.81 1155.81 0.00 
E076-10e 1135.46 1130.08 1117.77 -1.09 
E076-14s 1128.82 1122.8 1116.34 -0.58 
E101-08e 1428.8 1417.09 1391.74 -1.79 
El01-10¢ 1625.31 1605.11 1584.44 -1.29 
E101-14s 1550.85 1538.1 1512.92 -1.64 
Avg 953.79 941.59 938.14 -0.37 


表 4 货物 尺寸 信息 

Table 4 Size Information of Goods 
顾客 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
区 


1 13 16 16 了 Vb 6 L355 LY Tl le6 18 
MA 2 53 1 Oa Go .1 V5 Lb 1 i 3 1 9 1 


Var 和 汪汪 
Num 1 1 多 1 之 1 1 1 之 1 1 
1 1.2 12 人 91513113 1 23 
Ww 1 QB 15 TL LY 6 LY- 1 EB. 1 3 16 1 
h 1.2 Lo 2 O06 Oo 1 1 09 13 1 15 LY .1 i 
下 0 0 1 2 1 1 0 0 1 1 0 0 1 1 
wel 0.8 0.8 0.2 0.6 0.4 0.8 0.2 0.6 0.4 0.8 0.3 1.5 1.6 
Var 3 3 有 3 3 3 3 有 
Num i i 名 i i i 

1 1. L518 13 2 2.3 2 14 

MA i 3 15 O08 12 12 1 0.8 1 

h 1 1 刘 1 1 1 和 这 1 

和 i ' 0 0 lL 0 0 ' l 
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货车 采用 6.0x2.5x3.0 式 ， 空 载 质量 6t， 额 定 载重 9 t。 
货车 的 前 、 后 车 轴 位 置 的 最 大 承重 入 .N=4500kg 。 按 上 个 实 
验 ， 设 定 种 群 个 体 的 规模 50， 遗 传 迭 代 1 000 次 ， 交 叉 的 概 
率 0.8， 变 异 的 概率 0.05。 运 行 多 次 ， 取 最 好 结果 ， 见 图 4。 


a 


图 4 货车 行驶 路 径 
Fig.4 Vehicle traffic path 


图 4 可 知 ， 该 方案 有 四 条 送 货 路 径 ， 即 路 径 1: 


0-11-6-13-0;， 路 径 2: 0-1-12-15-9-0， 路 径 3: 0-3-8-5-10-0; 


本 路 径 4: 0-4-7-14-2-0。 货 车 行驶 的 总 路 程 为 311.3 km。 
区 在 同等 约束 下 ， 多 次 使 用 禁忌 算法 进行 实验 ， 取 最 好 结 
CN 果 。 松 弛 平衡 约束 , 进行 同样 的 实验 , 记录 数据 如 表 5 所 示 。 
© 表 中 dv 表示 结果 偏差 记 表示 最 好 运算 结果 江 表示 运行 时 间 。 
CC 表 5 与 禁忌 搜索 算法 比较 结果 

品 0) Table 5 Comparisons with tabu search algorithms 

LO 有 平衡 约束 无 平衡 约束 

[ew | P/km T/s P/km T/s 
©O) 多 阶段 算法 311.3 109.2 305.6 86.7 
OO 禁忌 算法 315.7 35.2 310.2 30.2 
QI dv/% -1.39 210.23 -1.48 187.09 
a 通过 对 表 5 横向 对 比 可 知 ， 平 衡 约束 不 仅仅 是 一 个 保证 
a 全 的 因素 ， 同 样 对 实验 结果 造成 影响 ， 是 一 个 不 可 忽略 的 
>< 约束 。 纵 向 进行 对 比 ， 多 阶段 算法 由 于 在 路 径 优 化 决策 时 ， 
(9 采用 了 回溯 算法 进行 深层 遍历 ， 导 致 实验 的 用 时 大 大 超过 了 
人 禁忌 搜索 算法 。 但 是 从 运算 的 结果 来 看 ， 多 阶段 算法 优 于 禁 
己 。 鼠 搜索 算法 。 表 明 本 文 提出 的 多 阶段 算法 在 解决 BL-VRP 上 
OO 有 效果 。 


4 ”结束 语 


本 文 首次 考虑 货物 的 重量 、 重 对 3L-VRP 进一步 扩 
展 ， 构 建 一 个 包含 纵向 轴 重 约束 、 横 向 重心 偏 移 约束 的 


BL-VRP 模型 。 针 对 BL-VRP 模型 ， 构 建 了 一 个 以 B-GA 算 
法 为 骨架 的 多 阶段 算法 。 通 过 两 组 实验 表明 ， 本 文 提 出 的 多 
阶段 算法 在 解决 3L-VRP 问题 上 不 逊 于 其 他 算法 ， 同 时 对 于 
BL-VRP 问题 也 表现 优秀 。 但 是 由 于 深层 搜索 的 使 用 ， 使 算 
法 的 时 效 性 较 差 ， 下 一 步 可 针对 此 处 进行 改进 。 
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